How to observe fluctuating temperature? 
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We provide arguments that event-by-event (EBE) analysis of multiparticle production 
data are ideal place to search for the possible fluctuation of temperature characterizing 
hadronizing source in thermodynamical approach. 

PACS numbers: 25.75.-q 24.60.-k 05.20.-y 05.70.Ln 

In this note we would like to advertise the feasibility and importance of the observations of 
fluctuating temperature in high energy multiparticle production processes. Let us first list our 
chain of reasoning: 

• ALICE will look (among other things) for the creation of Quark Gluon Plasma (QGP); 

• this calls for a thermodynamical description of (at least some) relevant observables and 
this in turn calls for the temperature T of the system under investigation as one of the 
most important quantities; 

• T can be deduced most directly (at least this is believed to be the case) by looking at px 



This is widely accepted approach, notwithstanding the fact that both details of what T 
really means or whether (|l]) is a proper form for /iy-dependence are still subject to hot debate 
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and modelling. Taking therefore (|I|) as our starting point we want to concentrate on question: 
is it possible that T is fluctuating quantity fl], ^ and if so, does it fluctuates only from event- 
to-event or also in a given event? 



We are not going to discuss here the problem of internal consistency (or inconsistency) of 
the notion of fluctuations of temperature in thermodynamics referring in this matter to [0, ^ 
What we want to do is to brings to ones attention the fact that event-by-event analysis allows 
(at least in principle) to detect fluctuations of temperature taking place in a given event. This 
is more than indirect measure of fluctuations of T proposed some time ago in or more direct 
fluctuations of T from event to event discussed in [|l|. 

To this end let us first remind our results presented in [0 where we have shown that 
fluctuation of the parameter in exponential distribution leads in a natural way (under some 
circumstances, of course, to be mentioned in a moment) to final distribution of the power-like 
form (know also as Levy distribution): 



exp 



l + (l-g)(^)-/iT 
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Here averaging (...) is performed over fluctuations of parameter (here our temperature) j; 
which take place around some mean value Tq and should follow gamma distribution (see 0] 
for details). The new parameter occuring here (identical to the so called entropic index or 
nonextensivity parameter in Tsallis statistics [§]) is tightly connected with the size of such 
fluctuations, namely 

q = I + u (3) 



where 



u = 



T / 

(4) 



It is worth to mention that distribution of the Levy type (^ has been observed already 
in inclusive processes 0. However, inclusive processes are not able to provide unambiguous 
answer what is the source of such behaviour. This can be done, such is our belief, only in the 
careful analysis of event-by-event data, especially those for heavy ion collisions. Two scenarios 
are possible here and should be subjected to experimental verification: 

(1) T is constant in each event but because of different initial conditions it fluctuates from 
event to event. In this case in each event one should find exponential dependence (|I]) with 
T = Tg^ent and possible departure from it will occur only after averaging over all events. 
It will reflect fluctuations originating in different initial conditions for each collision from 
which given event originates. This situation is illustrated in Fig. 1 where pt distributions 
for T = 200 MeV (black symbols) and T = 250 MeV (open symbols) are presented. All 
other details are the same as hsted below for Fig. 2. Such values of T correspond to 



2 



typical uncertainties in T expected at LHC due to different initial conditions. Notice that 
both curves presented here are straight lines. 



(2) T ffuctuates in each event around some value Tq. In this case one should observe departure 
from the exponential behaviour already on the single event level which should be fully 
given d^) with q > 1. It reflects situation when, due to some intrinsically dynamical 
reasons, different parts of a given event can have different temperatures [0]. Fig. 2 shows 
typical event of this type obtained in simulations performed for central Ph + Ph collisions 
taking place for beam energy equal i^team = 3 A-TeV in which density of particles in 
central region (defined by rapidity window — 1.5 < y < 1.5) is equal to ^ = 6000 (this 
is the usual value given by event generators like VENUS, SHAKER, HIJING). Black 
symbols represent exponential dependence obtained for T = 200 MeV (the same as in 
Fig. 1), open symbols show the power-like dependence as given by (0) with the same 
T and with q = 1.05 (notice that the corresponding curve bends slightly upward here). 
In this typical event we have ~ 18000 secondaries, i.e., practically the maximal possible 
number. Notice that points with highest pt correspond already to single particles. 

One should stress here the following important fact: our u = q — 1 has physical meaning of 
the total heat capacity C, because according to a basic relation of thermodynamic {P = -f) 

^ = ^ = - = '^-1- (5) 

Therefore measuring in addition to the temperature T also nonextensivity q describing its fluc- 
tuation (and, because of this, the total heat capacity C) could be of great practical importance 
for our understanding of dynamics of heavy ion collisions |l|, 0. In particular it should not only 
facilitate checking the commonly made assumption that an approximate thermodynamics state 
is obtained in a single collision but also, by knowing the heat capacity, provide considerable 
information about its thermodynamics (especially on the existence and type of the possible 
phase transitions ll|, @]) ||T0| . 
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Figure 1 
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Figure 2 
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